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ABSTRACT 

Quasi-periodic X-ray brightness oscillations (QPOs) with frequencies ^ 1 kHz have 
now been discovered in more than a dozen neutron stars in low-mass X-ray binary 
systems using the Rossi X-ray Timing Explorer. There is strong evidence that the 
frequencies of the kilohertz oscillations are the orbital frequencies of accreting gas in 
nearly circular orbits around these stars. Some stars that produce kilohertz QPOs may 
have spin frequencies ^400 Hz. For spin rates this high, first-order analytic treatments 
of the effects of the star's rotation on its structure and the spacetime are inaccurate. 
Here we use the results of a large number of fully relativistic, self-consistent numerical 
calculations of the stellar structure of rapidly rotating neutron stars and the interior and 
exterior spacetime to investigate the constraints on the properties of such stars that can 
be derived if stable circular orbits of various frequencies are observed. We have computed 
the equatorial radius of the star, the radius of the innermost stable circular orbit, and 
the frequency of the highest-frequency stable circular orbit as functions of the stellar spin 
rate, for spin rates up to the maximum possible and for several illustrative equations 
of state. Our calculations show that the upper bounds on the stiffness of neutron star 
matter implied by a given orbital frequency are typically significantly stricter for stars 
with spin frequencies ^ 400 Hz than for slowly rotating stars. 

Subject headings: accretion, accretion disks — dense matter — equation of state — 
gravitation — relativity — stars: neutron 



1. INTRODUCTION 

The successful launch of the Rossi X-ray Timing 
Explorer (RXTE) has made it possible to investigate, 
for the first time, the X-ray variability of neutron stars 
and black holes at frequencies > 300 Hz. One of the 
most important discoveries made with RXTE is that 
many neutron stars in low-mass X-ray binaries pro- 
duce high-frequency brightness oscillations with fre- 
quencies in the range ~300-1200 Hz (see van der Klis 
1998). High-frequency oscillations are observed both 
during type I (thermonuclear) X-ray bursts and in the 
persistent X-ray emission. The discovery of these os- 
cillations has made possible derivation of interesting 
upper bounds on the masses and radii of these neu- 
tron stars and significant new constraints on the equa- 
tion of state of neutron star matter (Miller, Lamb, & 
Psaltis 1998a, 1998b; Lamb, Miller, & Psahis 1998; 
Miller & Lamb 1998; Strohmayer et al. 1998). 

Only a single oscillation has been observed from 
each source during a type I X-ray burst, and the os- 
cillations in the tails of bursts appear to be highly 
coherent (see, e.g., Smith, Morgan, & Bradt 1997), 
with frequencies that are always the same in a given 
source (see, e.g., Strohmayer 1997). The burst oscil- 
lations are thought to be caused by the existence of 
one or two brighter regions on the stellar surface that 
produce oscillations at the stellar spin frequency or 
its first overtone as the star spins (see Strohmayer et 
al. 1997b for compelling evidence in favor of this in- 
terpretation) . The frequency of the burst oscillations 
ranges from ~330 Hz in 4U 1702-42 (Swank 1997) 
to 589 Hz in an unidentified source in the direction of 
the Galactic center (Strohmayer et al. 1997a). 

The kilohertz quasi-periodic oscillations (QPOs) 
observed in the persistent emission have high ampli- 
tudes and relatively high coherences. The frequencies 
of the two QPOs often observed simultaneously in a 
given source have a frequency separation that is al- 
most constant in all sources (sec Wijnands & van der 
Klis 1997; Psaltis et al. 1998a, 1998b), although the 
frequencies of the QPOs themselves vary by hundreds 
of Hertz (see Wijnands ct al. 1998 and van der Klis 
1995). The separation frequencies of the two kilohertz 
QPOs seen in 4U 1728-34 (Strohmayer et al. 1996, 
1997b) and 4U 1702-42 (Swank 1997) arc consistent 
with the frequencies of their burst oscillations. The 
separation frequencies in 4U 1636—536 (Wijnands et 
al. 1997; Zhang et al. 1997) and KS 1731-260 (Smith 
et al. 1997; Wijnands & van der Klis 1997) are con- 



sistent with half the frequencies of their burst oscilla- 
tions. 

The presence of only two simultaneous kilohertz 
QPOs in a given source, the approximately constant 
frequency separation Ap between them, and the con- 
sistency of Al/ with the stellar spin frequency inferred 
from burst oscillations is strong evidence that the stel- 
lar spin is generating the frequency difference, that 
only one sideband of the primary QPO frequency is 
being generated, and that one of the two QPOs is 
therefore caused by the beat of the spin frequency 
with the other frequency (see Lamb et al. 1998; Miller 
et al. 1998a). This implies that in addition to the 
spin frequency there is only one other primary fre- 
quency and that this frequency is a rotational fre- 
quency such as an orbital frequency. This excludes 
neutron star surface and photon bubble oscillations 
as explanations for the primary kilohertz QPO fre- 
quency and makes disk oscillations an improbable ex- 
planation (see van der Klis 1998; Lamb ct al. 1998). 

In all sources the frequencies of the kilohertz QPOs 
fall within the expected range of orbital frequencies 
near a neutron star and can vary by several hundred 
Hertz in less than a few hundred seconds while re- 
maining highly coherent [yj Av ^ 100). This is strong 
further evidence against disk oscillations and in favor 
of orbital motion of inhomogeneities in the accretion 
disk as the cause of the primary kilohertz QPO (see 
Lamb et al. 1998). The accreting gas exerts a strong 
torque on the neutron star and hence the star is ex- 
pected to be spinning in the same sense as the orbital 
motion of the accreting gas. The beat frequency must 
therefore be the lower of the two kilohertz QPO fre- 
quencies whereas the orbital frequency is the higher. 

There are two candidates for the orbital frequency: 
the frequency at the radius where the accreting gas 
first couples strongly to magnetic field of the neutron 
star and the frequency at the sonic radius where ra- 
diation forces or general relativistic effects cause the 
radial motion of the gas to increase sharply and be- 
come supersonic. The general properties of the kilo- 
hertz QPO sources and the specific properties of the 
kilohertz QPOs themselves strongly indicate that the 
relevant frequency is the orbital frequency at the sonic 
point (Miller et al. 1998a; see also Lamb et al. 1998). 
In either case, the frequency t'Qpo2 of the higher- 
frequency kilohertz QPO is the orbital frequency of 
gas in a nearly circular orbit around the neutron star, 
whereas the frequency fQpoi of the lower-frequency 
QPO is the beat of the neutron star spin frequency 
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with an orbital frequency near ;^qpo2- 

In the following discussion we shall describe a cir- 
cular orbit as stable or unstable according to its prop- 
erties as determined by solving the geodesic equation 
for a test particle moving in that orbit in the space- 
time of interest. However, it is important to bear 
in mind that there are no closed, circular orbits in 
the vicinity of an accreting neutron star, because the 
motion of gas near such a star is affected not only by 
the curvature of spacetime but also by radiation, mag- 
netic, and viscous forces, which cause the gas to spiral 
inward even at radii where, in their absence, closed, 
stable orbits would be possible (Miller & Lamb 1993, 
1996). However, the distinction between stable circu- 
lar orbits (SCOs) and unstable circular orbits is still 
relevant. In particular, the innermost stable circular 
orbit (ISCO) is still physically significant when the ef- 
fects of radial gas pressure forces on the ISCO can be 
neglected (which should be valid when the luminosity 
of the source is much less than the Eddington critical 
luminosity), because under these conditions gas inside 
the ISCO spirals inward so quickly that it cannot pro- 
duce a wavetrain with the coherence observed for the 
kilohertz QPOs, regardless of whether it is acted on 
by radiation, magnetic, and viscous forces (Miller et 
ai. 1998a). 

Identification of the higher-frequency kilohertz QPO 
with the frequency of an SCO has made it possi- 
ble to derive upper bounds on the masses and radii 
of the neutron stars in the kilohertz QPO systems 
(Miller et al. 1998a, 1998b). These bounds follow 
from the requirement that the radius i?orb of the or- 
bit be larger than the radius i?nis of the ISCO as well 
as larger than the equatorial radius i?eq of the star (if 
Rcq < i?orb < ^ms, orbits with the required frequency 
exist but are unstable). 

For nonrotating stars, observation of a given or- 
bital frequency can be used to derive upper bounds on 
the mass and radius that are independent of the equa- 
tion of state assumed (Miller et al. 1998a). For rotat- 
ing stars, the situation is more complicated. In gen- 
eral, both the structure of the star and the spacetime 
are affected by the star's rotation, and there are no 
general analytical expressions for the relevant quanti- 
ties. However, the exterior spacetime of a slowly and 
uniformly rotating fluid star is unique to first order in 
the dimensionless angular momentum j = cJ/GM'^, 
where J and M are the star's angular momentum 
and gravitational mass, and can be expressed ana- 
lytically to this order (Hartle & Thorne 1968). The 



leading c;orrections to the expressions for the orbital 
frequency and the radius of the ISCO are first-order 
in j. Using these expressions, one can demonstrate 
that observation of a given orbital frequency also im- 
plies upper bounds on the mass and radius of a slowly 
rotating star (Miller et al. 1998a). For a given stel- 
lar spin frequency, these upper bounds depend on the 
moment of inertia and hence on the equation of state 
assumed. 

Many of the kilohertz QPO sources appear to 
have spin frequencies ~250-350 Hz (see Miller et al. 
1998a). Examples of such sources include 4U 0614-1-091, 
4U 1608-52, 4U 1820-30, Cyg X-2, Sco X-1, GX 5-1, 
and GX 17-1-2, all of which have kilohertz QPO sepa- 
ration frequencies in this range, as well as 4U 1728—34 
and 4U 1702-42, which not only have kilohertz QPO 
separation frequencies in this range but also have 
burst oscillation frequencies that are consistent with 
these separation frequencies. The ~520 and ~580 Hz 
frequencies of the burst oscillations seen in KS 1731—260 
and 4U 1636—536 are thought to be twice their spin 
frequencies, although this is not certain. For spin fre- 
quencies in this range, j is ~0.1-0.3, depending on 
the assumed equation of state and the mass of the 
star, and hence an analysis that is first-order in j is 
quite accurate for such stars. Such an analysis shows 
that spin rates ~300 Hz can increase the upper bound 
on the stellar mass by as much as ~10-20% but typi- 
cally increase the upper bound on the radius by only 
~2-5% (Miller et al. 1998a, 1998b). 

On the other hand, some neutron stars that show 
kilohertz QPOs may turn out to have spin frequen- 
cies > 400 Hz. For example, oscillations with frequen- 
cies ~550 and ~590 Hz have been seen during X-ray 
bursts from, respectively, Aql X-1 and the unknown 
source in the direction of the galactic center, possibly 
indicating that these neutron stars have spin frequen- 
cies this high (Miller et al. 1998a). The recent discov- 
ery using RXTE that the source SAX J1808.4-3658 
has a coherent 401 Hz oscillation indicates that this 
accreting neutron star is spinning rapidly (Wijnands 
& van der Klis 1998a, 1998b). For stars spinning this 
fast, the effect of the star's spin on its equilibrium 
structure (which is second-order) and on the space- 
time can be substantial. In order to obtain accurate 
results for such high spin rates, the equilibrium stel- 
lar structure and the interior and exterior spacetime 
must be computed self-consistently, which can only 
be done numerically. 

Here we use the results of a large number of fully 
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rclativistic, sclf-consistciit numerical calculations of 
the structure of rapidly rotating neutron stars and the 
interior and exterior spacetime to investigate the con- 
straints on the properties of such stars that can be de- 
rived if SCOs of various frequencies are observed. We 
have computed the equatorial radius of the star, the 
radius of the ISCO, and the frequency of the highest- 
frequency SCO as functions of the stellar spin rate 
and gravitational mass, for spin rates up to the maxi- 
mum possible and for several illustrative equations of 
state. Comparison of these results with the highest 
observed kilohertz QPO frequency in a given source 
can be used to derive bounds on the mass and ra- 
dius of the neutron star in that source, for a given 
equation of state. We also report the frequency of 
the highest-frequency SCO as a function of the stel- 
lar spin rate, for stars of any mass constructed using 
a given equation of state. These curves can be used 
to check whether a particular equation of state is con- 
sistent with the frequency of a given kilohertz QPO. 

Our calculations show that the upper bounds on 
the stiffness of neutron star matter implied by a given 
orbital frequency are typically significantly stricter for 
stars with spin frequencies > 400 Hz than for slowly 
rotating stars. 

In § 2 we describe our assumptions and calcula- 
tional methods. In § 3 we present our results and dis- 
cuss the implications for constraining the properties 
of neutron star matter. Our conclusions are summa- 
rized in § 4. 

2. ASSUMPTIONS AND METHOD 

In deriving bounds on the masses and radii of the 
neutron stars with kilohertz QPOs, we assume that 
the higher-frequency of the two simultaneous kilo- 
hertz QPOs is the frequency of a stable circular orbit 
around the neutron star, for the reasons discussed in 
§ 1. Hence the orbital radius i?orb that corresponds 
to the QPO frequency must be larger than both the 
equatorial radius i?eq of the neutron star and the ra- 
dius i?ms of the ISCO (Miller et al. 1998a). 

We compute the equilibrium stellar structure and 
the interior and exterior spacetime using the numer- 
ical code described in Cook, Shapiro, & Teukolsky 
(1992, 1994a, 1994b). This code solves the full gen- 
eral relativistic equation of hydrostatic equilibrium 
for a star with a given spin rate using a variation of 
the metric potential method of Komatsu, Eriguchi, & 
Hachisu (1989a, 1989b). It gives accurate solutions 



even for stars that are spinning very rapidly. 

In any stationary, axisymmetric spacetime, the or- 
bital frequency at a given coordinate radius, as mea- 
sured at infinity, is = d(j)/dt, where (j) and t are, 
respectively, the global azimuthal and time coordi- 
nates based on the spacelike and timelike Killing vec- 
tor fields of the spacetime. The time interval required 
for one orbit of an element of gas is the same ev- 
erywhere, as measured in the global time coordinate. 
Given the metric of the exterior spacetime, the or- 
bital frequency at a given radius is the solution of the 
equation (see Lightman et al. 1973, p. 469) 

g<l,<t,,r^'^ + 2gtfj,^r^ + gtt,r = , (1) 

where g^^, gt0, and gu are the metric components 
indicated and commas denote partial derivatives. 

2.1. Masses and Equations of State 

We have explored the constraints implied by ob- 
servation of an SCO of given frequency for a variety 
of neutron star equations of state. These restrictions 
are most significant if the equation of state is hard 
rather than soft. Hence, in this report we present re- 
sults for four relatively hard equations of state. For 
completeness, we consider both baryonic masses that 
are stable for nonrotating stars (the so-called "nor- 
mal" sequences of Cook et al. [1994b]) and the higher 
baryonic masses that are stable only for rotating stars 
(the "supramassive" sequences of Cook et al. [1994b]). 
Whether the supramassive sequences are accessible 
depends on how the specific angular momentum of 
the accreting gas evolves with time. 

In order to facilitate comparisons with previous 
studies of neutron star properties (see, e.g.. Pethick & 
Ravenhall 1995), we consider the Friedman-Pandhari- 
pande-Skyrme (FPS) equation of state (Friedman & 
Pandharipande 1981; Lorenz, Ravenhall, & Pethick 
1993). The FPS equation of state is based on the 
Urbana v\i two-nucleon potential plus the density- 
dependent three- nucleon interaction model of Lagaris 
& Pandharipande (1981) and gives a maximum gravi- 
tational mass for a nonrotating star of about 1.8 Mq, 
compared with a maximum mass of 2.12 Mq for a ro- 
tating star. The maximum spin frequency for stars in 
the normal sequence is 1411 Hz, and the maximum 
spin frequency for stars in the supramassive sequence 
is 1878 Hz. 

As an example of later realistic equations of state, 
we consider the UU equation of state (Wiringa, Fiks, 
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& Fabrocini 1988), which is based on the Urbana 
Vi4 two-nucleon potential plus the Urbana VII three- 
nucleon potential (Schiavilla, Pandharipande, & Wiringa 
1986) and gives a maximum mass for a nonrotating 
star of about 2.2Mq. Although it is based on older 
scattering data, the UU equation of state is similar to 
the recent A18 + UIX' + 5ub equation of state (Akmal, 
Pandharipande, & Ravenhall 1998), which is based on 
the modern Argonne vis two-nucleon potential and 
the Urbana IX three-nucleon potential and takes into 
account the nonzero momentum of the interacting 
nucleons (see Pandharipande, Akmal, & Ravenhall 
1998). Like the A18 + UIX' + Svh equation of state, 
the UU equation of state gives a maximum mass of 
about 2.2 Mq for a nonrotating neutron star. The 
maximum mass for a rotating neutron star is 2.61 Mq, 
and the maximum rotation frequencies for the normal 
and supramassive sequences are, respectively, 1561 Hz 
and 1989 Hz. 

In order to illustrate the generic effects of signifi- 
cant softening of a stiff equation of state at a critical 
density, we consider the tensor interaction (TI) equa- 
tion of state of Pandharipande & Smith (1975a; M in 
the Arnctt & Bowers [1977] survey). Although the TI 
equation of state is itself no longer of interest to nu- 
clear physicists, this equation of state demonstrates 
the effects of a very strong, first-order phase transi- 
tion, such as may occur at the transition from nu- 
cleon matter to quark matter (see Glendenning 1992; 
Heiselberg, Pethick, & Staubo 1993; Pandharipande 
et al. 1998). The maximum mass of a nonrotating star 
constructed using the TI equation of state is 1.8 Mq, 
and the maximum mass of a rotating star is 2.1 Mq. 
The maximum rotation frequency for the normal se- 
quence is 707 Hz, and for the supramassive sequence 
is 1229 Hz. 

Finally, as an example of the relatively stiff equa- 
tions of state often given by mean field theories, we 
consider the mean-field equation of state of Pandhari- 
pande & Smith (1975b; L in the Arnett & Bowers 
[1977] survey). The maximum mass of a nonrotating 
star constructed using this equation of state is 2.7M0, 
compared with 3.27Mq for a rotating star. The max- 
imum rotation frequency is 1031 Hz for the normal 
sequence and 1321 Hz for the supramassive sequence. 

2.2. First-Order Expressions 

In § 3 it will be instructive to compare the be- 
havior of the orbital frequencies and radii computed 
using our numerical models of rapidly rotating stars 



with the behavior given by the analytical expressions 
valid for slowly rotating stars. As noted in § 1, the 
spacetime around a rotating fluid star is unique to 
first order in the dimensionless angular momentum j. 
To this order in j, the frequency of a prograde orbit 
at circumferential radius r around a star with gravi- 
tational mass M is (see Lightman et al. 1973, p. 469; 
Miller et al. 1998a) 



n = 



l-j{M/rf^^] (M/r^) 



and the circumferential radius of the ISCO is 



R^s{M,j)^6M[l~j{2/3f/^] 



(2) 



(3) 



in units in which G = c = 1. In the present work 
we always quote circumferential radii, in contrast 
to Miller et al. (1998a), where we quoted Boyer- 
Lindquist radii. The two radii are identical to first 
order in j, but to higher orders in j the circumferen- 
tial radius is the physically meaningful radius, which 
is the reason we use it here. Combining equations (^ 
and one can show that to first order in j, the fre- 
quency i^K,ms of the innermost stable prograde orbit 
is (see Kluzniak, Michelson, & Wagoner 1990; Miller 
et al. 1998a) 



2210 (1 -I- 0.75j)(Mo/M) Hz 



(4) 



Thus, for slowly rotating stars the frequency of the 
ISCO increases linearly with the star's spin rate. 

Using expressions (^~(|^), one can show (Miller et 
al. 1998a) that the mass and radius of a slowly rotat- 
ing star are bounded above by 



M^ax« [l+0.75j(;^spin)]M° 



and 



R„ 



[1 + 0.20j>,pin)]<, 



(5) 



(6) 



Here j(i'spin) is the value of j for the observed stel- 
lar spin rate at the maximum allowed mass for the 
equation of state being considered and 



and 



AC = 2.2 (1.0 kHz/z.*po2) Mg 
i?^,, = 19.5 (1.0 kRz/^*^po2) km 



(7) 



(8) 



I non- 



^-max i^.^ v^-" / -i^fu^J 

are the upper bounds on the mass and radius of a ^^^^^ 
rotating star in terms of i^Qpo2i the highest observed 
frequency of the higher-frequency kilohertz QPO. The 
precise upper bounds on the mass and radius depend 
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on the equation of state and can be determined by 
searching a grid of neutron star models for the one 
that gives the maximum allowed mass. Equations (^) 
and (^) show that the bounds are always greater for 
a slowly rotating star than for a nonrotating star, re- 
gardless of the equation of state assumed. 

No expressions similar to equations (||)-(||) are 
available for rapidly rotating stars. 

3. RESULTS AND DISCUSSION 

We first show how the radius of the ISCO and the 
equatorial radius vary with stellar spin rate for stars 
constructed using the FPS equation of state. The be- 
havior of these radii makes clear why the frequency 
of the highest-frequency SCO around a star of given 
mass generally first increases as the star is spun up 
and then decreases. Considering this behavior for sta- 
ble stars with different masses makes the behavior of 
the maximum frequency of an SCO for stars of any 
mass understandable. 

Next, we present mass-radius relations for stars 
with a wide range of spin rates, constructed using the 
FPS and UU equations of state. We then show how 
to derive limits on the mass and radius of a rapidly 
rotating star from the frequency of an SCO around it 
and discuss the constraints on the equation of state 
of neutron star matter implied by such constraints. 

3.1. Radii and Orbital Frequencies 

Figures la and lb show how the circumferential 
radius of the ISCO and the circumferential radius of 
the stellar equator vary with I'spim the stellar spin fre- 
quency measured at infinity, for stars constructed us- 
ing the FPS equation of state. These stars have con- 
stant baryonic masses equal to those of non-rotating 
stars with gravitational masses of IAMq and 1.6Mq. 
Because the gravitational mass increases only slightly 
with increasing spin frequency, the curves for stars 
of constant gravitational mass are almost identical 
(the largest stable equatorial radii are very slightly 
smaller). As expected, the dimensionless angular 
momentum j increases linearly with spin rate for 
slowly rotating stars but more steeply for rapidly 
rotating stars: for the 1.4 A/q model, j = 0.23 at 
I'spin = 0.5 kHz and 0.52 at 1.0 kHz; for the 1.6 Mg 
models j = 0.20 at 0.5 kHz and 0.43 at 1.0 kHz; for 
the 1.8 Mq models j = 0.16 at 0.5 kHz and 0.36 at 
1.0 kHz. Figure Ic shows how the frequency of the 
highest-frequency SCO varies with stellar spin rate of 



1.4 Mq, 1.6 Mq, and 1.8 Af© stars. 

The circumferential radius i?nis of the ISCO de- 
creases linearly with spin rate for slowly rotating 
stars, in agreement with the first-order expression (^, 
but decreases more slowly as the spin rate increases. 
For the IAMq star, the deviation from equation (||) 
is significant at ^spin ~ 300 Hz (see Fig. la). For the 
1.6 Mq star, the deviation is significant at 500 Hz and 
at about 1065 Hz i?nis reaches a minimum and then 
increases with increasing spin rate, until an ISCO no 
longer exists (see Fig. lb). In contrast, the circumfer- 
ential equatorial radius of the stellar models increases 
quadratically with the spin rate from i^spin = Oj ex- 
ceeding i?ms at about 580 Hz for the 1.4 A/q star 
and at about 1220 Hz for the 1.6 ATq star. For spin 
rates above these critical rates, all circular orbits with 
radii larger than the star's equatorial radius are sta- 
ble: there is no ISCO. 

The Kerr spacetime can be expressed analytically 
and is therefore sometimes used as a convenient ap- 
proximation to the exterior spacetime of a spinning 
neutron star. For this reason, in Figures la and lb we 
compare i?nis(Kerr), the circumferential radius of the 
ISCO in a Kerr spacetime with the same gravitational 
mass and angular momentum as the stellar models, 
with the actual radius i?ins of the ISCO. Unlike the 
actual radius, i?nis(Kerr) decreases monotonically and 
nearly linearly with spin rate even at high spin fre- 
quencies. Indeed, at high spin rates i?nis(Kerr) de- 
creases faster than linearly with increasing spin rate, 
and hence the exact Kerr expression for Rj^s is a worse 
approximation than the first-order approximation (3). 

i?nis(Kerr) is significantly smaller than R^s at high 
spin rates. As a result, when R^s is larger than 
the equatorial radius of the star, the frequency of 
the highest-frequency SCO is significantly lower than 
one would estimate using the Kerr spacetime, and the 
constraints on the mass and radius of the star are cor- 
respondingly tighter. For both the IAMq and 1.6AfQ 
stars, the critical spin rate at which the ISCO disap- 
pears in the Kerr approximation is about 23% smaller 
than in the actual spacetime. Thus, for stellar spin 
rates > 400 Hz, the exterior spacetime of a spinning 
black hole is generally an inaccurate approximation 
to the exterior spacetime of a neutron star. 

Figure Ic shows why the constraints on the mass 
and radius of a slowly rotating star implied by a given 
SCO frequency are generally looser for a slowly rotat- 
ing star than for a nonrotating star of the same mass, 
whereas the constraints on a rapidly rotating star 
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Fig. 1. — Typical variations of important radii and 
frequencies with stellar spin rate, (a) Circumferential 
radii of the innermost stable circular orbit (dashed 
line) and the stellar equator (solid line) as a function 
of spin rate, for a \AMq star. Also shown for com- 
parison is the circumferential radius of the innermost 
stable circular orbit in a Kerr spacetime with the same 
gravitational mass and angular momentum (dotted 
line), (b) Same radii as in (a), but for a I.6M0 star, 
(c) Frequency of the highest-frequency stable circu- 
lar orbit as a function of stellar spin rate, for 1.4 M© 
(dotted hue), 1.6 Mq (dashed Hue), and 1.8 M© (solid 
dashed line) stars. All stellar models were constructed 
using the FPS equation of state. 



are usually much tighter. For slowly rotating stars 
with gravitational masses of 1.4, 1.6, and 1.8 M© con- 
structed using the FPS equation of state, the equato- 
rial radius of the star is smaller than the radius of the 
ISCO. Hence, at low spin rates the highest-frequency 
SCO is the ISCO, which at these spin rates shrinks 
linearly as the spin rate increases (see eq. Q), caus- 
ing the frequency of the highest-frequency SCO to 
increase linearly with the star's spin rate. However, 
at a certain critical spin rate the equatorial radius of 
the star becomes larger than the radius of the ISCO 
for a star with the given gravitational mass and angu- 
lar momentum. For spin rates above this critical spin 
rate, the highest-frequency SCO is the orbit that just 
skims the stellar surface. At high spin rates, the equa- 
torial radius of the star increases rapidly with increas- 
ing spin rate, causing the frequency of the highest- 
frequency SCO to decrease rapidly. 

For the M — 1.4 M© star, the frequency of the 
highest-frequency SCO is maximized at the spin fre- 
quency for which the radius of the ISCO is equal to 
the radius of the stellar equator (see Fig. la). How- 
ever, the highest-frequency SCO for a star of given 
mass but any possible spin rate is not necessarily the 
ISCO with radius equal to the equatorial radius of the 
star. This is illustrated by the M = 1.6 Af© star. For 
this star, the frequency of the highest-frequency SCO 
has its maximum at the spin frequency at which the 
radius of the ISCO is a minimum (see Fig. lb). At 
this frequency the radius of the ISCO is larger than 
the equatorial radius of the star. 

3.2. Maximum SCO Frequency 

The maximum frequency of the highest-frequency 
SCO for stable stars of any mass can be determined 
by constructing curves like those shown in Figure Ic, 
for a dense sequence of stellar masses. The curve of 
maximum frequency as a function of spin rate is then 
the upper envelope of these curves. Figure 2 shows 
curves of maximum SCO frequency as a function of 
stellar spin rate, for stable stars of any mass con- 
structed with the four indicated equations of state. If 
the measured spin frequency of the star and the fre- 
quency of a nearly circular orbit correspond to a point 
that lies above the curve for a given equation of state, 
that equation of state is excluded for all neutron stars. 

As Figure 2 shows, the maximum frequency of the 
highest-frequency SCO /or stars of any mass typically 
decreases with increasing spin rate, even though the 
frequency of the highest-frequency SCO for a fixed 
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Fig. 2. — Maximum frequency of a stable circular 
orbit as a function of spin rate for stable stars of 
any mass, for the four indicated equations of state, 
which are discussed in the text. The curve labeled 
"Mnormai" iucludes Only normal sequences, whereas 
the other curves include both normal and supramas- 
sive sequences (see § 3.1 for a discussion). 

gravitational mass increases with increasing spin rate 
over a wide range of spin rates. The reason is that 
at low spin rates the mass that gives the highest- 
frequency SCO is the mass at which the ISCO co- 
incides with the stellar equator. This mass increases 
with increasing spin rate, causing the maximum fre- 
quency to decrease. 

Stars constructed with the TI equation of state M 
are an exception to this rule. This equation of state 
is extremely stiff at the low densities characteristic 
of the centers of lower-mass stars, but becomes very 
soft at the density reached at the center of a nonro- 
tating M k, 1.75 star, because a pion condensate 
forms at this density. For this equation of state, the 
maximum orbital frequency occurs for a stellar mass 
near the maximum mass allowed by the star's spin 
frequency (e.g., \.%Mq for a slowly rotating star). 

Consider first the normal sequences. At low spin 
frequencies, the surface of the 1.8 star is well in- 
side the ISCO, the highest-frequency SCO is there- 
fore the ISCO, and the frequency of the maximum- 



frequency SCO therefore increases with increasing 
spin rate. However, at about 500 Hz the maximum 
frequency stops increasing and then plummets, as 
shown by the curve labeled "Mnormar in Figure 2. 
The reason is a general rclativistic effect first pointed 
out by Cook et al. (1994b): as the angular momen- 
tum of a star increases, the spin frequency first in- 
creases, then decreases, and finally increases again, 
producing a local maximum in the spin frequency vs. 
angular momentum relation. For normal sequences 
using equation of state M, this local maximum oc- 
curs at a spin frequency slightly greater than 500 Hz. 
Hence, in order to have an observed spin frequency 
higher than this, the star must have a much higher 
angular momentum, but a star with this much angu- 
lar momentum has a much larger equatorial radius, 
larger than the ISCO for a star of its mass and an- 
gular momentum, so the highest-frequency SCO is at 
the stellar surface and has a smaller frequency. 

Consider now the supramassive sequences. For a 
fixed baryonic mass, there is a local maximum in the 
spin frequency vs. angular momentum relation, just 
as for the normal sequences. A curve of I'sco.max vs. 
J'spin constructed for a supramassive sequence with a 
fixed baryonic mass would therefore look similar to 
the "Mnormai" curve in Figure 2, except that it would 
start at a positive spin frequency and would plum- 
met at a higher spin frequency. The local maximum 
of the spin frequency increases with increasing bary- 
onic mass. The envelope of these curves produces 
the curve labeled "M". The rapid downturn of this 
curve at a spin frequency of approximately 1100 Hz 
occurs because above this frequency there is no ISCO, 
and consequently the highest-frequency orbit is the 
one that skims the surface. For Vspi-a > HOO Hz, 
the equatorial radius increases rapidly with increasing 
spin frequency, and therefore the maximum frequency 
of a circular orbit decreases rapidly. 

Comparison of the two curves in Figure 2 for equa- 
tion of state M demonstrates the potential impor- 
tance of the supramassive scqiiences. For example, if 
only normal sequences are physically accessible (e.g., 
if the gas accretes with low specific angular momen- 
tum), then observation of an 1100 Hz SCO from a star 
with a spin frequency in excess of 550 Hz would rule 
out equation of state M, whereas equation of state 
M would still be viable if supramassive sequences are 
accessible. 



8 




Fig. 3. — Constraints on neutron star masses and radii imposed by stellar stability and observation of a stable 

circular orbit, (a) Solid lines show the niass-radius relations for FPS stars with the spin frequencies indicated 
(in kHz). The high- mass end of each mass-radius relation shown is the radial instability limit (short-dashed line) 
whereas the low-mass end is the mass shedding limit (dotted line). The dash-dotted line shows the lowest stellar 
mass (largest stellar radius) consistent with the requirement that the radius of the star be smaller than the radius of 
a 1200 Hz orbit. For FPS stars and this orbital frequency, the requirement that the radius of the orbit also exceed 
the radius of the ISCO does not constrain the mass or radius of the star. The bold portion of each mass-radius 
curve highlights the region allowed by both the physical limits of the equation of state and by the observation of an 
SCO with the listed frequency, (b) Same as in (a) but for UU stars and an orbital frequency of 1400 Hz. Although 
the dash-dotted line curves strongly to the left at high stellar spin rates, it always shows the smallest stellar mass 
(largest stellar radius) consistent with the requirement that the radius of the star be smaller than the radius of a 
1400 Hz orbit. The line of long dashes shows the largest stellar mass (smallest stellar radius) consistent with the 
requirement that the radius of the orbit exceed the radius of the ISCO. 



3.3. Mass cind Radius Bounds 

As explained in § 1, observation of a given SCO fre- 
quency around a nonrotating star allows one to derive 
upper bounds on the mass and radius that are inde- 
pendent of the equation of state, whereas for a rotat- 
ing star one must consider a specific equation of state 
in order to derive bounds on the mass and radius. 

Given an equation of state and a stellar spin rate, 
the mass of the star must be between the mass- 
shedding limit and the radial instability limit, regard- 
less of the values of any orbital frequencies. The ra- 
dius of the star is bounded by the extreme values of 
the radii given by the mass-radius relation over this 



mass interval. Observation of an SCO with a cer- 
tain frequency may allow one to restrict further the 
allowed mass and radius intervals, depending on the 
equation of state and the frequency of the SCO. The 
possible further restrictions are of two types: a lower 
bound on the mass, imposed by the requirement that 
the radius of the orbit be greater than the radius of 
the star, and an upper bound on the mass, imposed by 
the requirement that the radius of the orbit be equal 
to or greater than the radius of the ISCO. If either of 
these bounds restrict further the allowed mass range, 
the radius of the star is bounded by the extreme val- 
ues of the radii given by the mass-radius relation over 
this reduced mass interval. 
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Figure 3 shows the constraints on neutron star 
masses and radii imposed by stellar stability and ob- 
servation of a stable circular orbit, for stars con- 
structed using the FPS and UU equations of state (see 
§ 2.2). The mass-radius relations (curves of gravita- 
tional mass versus equatorial circumferential radius) 
shown in this figure were constructed by generating 
several sequences of stellar models. Each sequence 
consisted of stellar models with the same baryon num- 
ber but a range of spin frequencies. The grid of mod- 
els constructed in this way was then used to generate 
the mass-radius relations shown in Figure 3. These 
relations arc tabulated in Tables 1 and 2, where for 
each spin frequency the boundary between normal 
and supramassive stars is indicated by the dashed 
lines. The relations shown for rapidly spinning stars 
are much flatter than the usual mass-radius relations 
for nonrotating stars. 

The high- mass end of each constant spin-frequency 
sequence shown in Figure 3 is the gravitational mass 
above which the star is unstable to a radial instability. 
This mass limit is indicated by the short-dashed line 
in Figures 3a and 3b. The maximum gravitational 
mass increases with increasing spin rate, both because 
the gravitational mass corresponding to a given bary- 
onic mass increases with increasing spin and because 
the maximum stable baryonic mass increases with in- 
creasing spin. The low-mass end of each mass-radius 
relation shown is the gravitational mass below which 
the star is subject to mass-shedding at the equator. 
This mass limit is indicated by the dotted lines in 
Figures 3a and 3b. 

Consider now the possible further restrictions on 
the allowed mass and radius intervals imposed by ob- 
servation of an SCO with a high frequency. The re- 
quirement that the radius of an orbit be larger than 
the equatorial radius of the star places a lower bound 
on the mass of the star given the frequency of an SCO. 
If, for that star's spin rate, there are stable stars with 
masses smaller than this lower bound, then the ob- 
servation of an SCO raises the lower bound on the 
mass of a star. Whether observation of an SCO with 
a given frequency raises the lower bound on the stel- 
lar mass depends on the equation of state and the 
star's spin rate as well as the frequency of the SCO. 
For example, for the FPS equation of state and a spin 
rate 600 Hz, observation of a 1200 Hz SCO imposes 
a lower bound on the mass of 0.8 M©, whereas the 
mass-shedding limit is 0.56 Mq (see Fig. 3a). Hence, 
in this case the limit imposed by the SCO is stricter. 



On the other hand, if z^spin > 1200 Hz , there is always 
an SCO with a frequency of at least 1200 Hz around 
any FPS star that is stable against mass-shedding and 
hence the SCO observation does not further restrict 
the allowed mass interval. In contrast, the lower mass 
limit imposed by observation of an SCO with a fre- 
quency of 1400 Hz is always stricter than the mass 
shedding limit for a UU star, regardless of its spin 
rate. 

Observation of an SCO lowers the upper bound on 
the mass of a star if, for that star's spin rate, there 
are stable stars with ISCOs with radii larger than that 
permitted by the requirement that the radius of the 
orbit be larger than the radius of the ISCO. Again, 
whether observation of an SCO with a given frequency 
lowers the upper bomid on the stellar mass depends 
on the equation of state and the star's spin rate as 
well as the frequency of the SCO. For example, for the 
FPS equation of state observation of a 1200 Hz SCO 
frequency would not lower the upper bound on the 
stellar mass imposed by the radial instability limit, 
regardless of the star's spin rate. In contrast, for the 
UU equation of state observation of a 1400 Hz SCO 
frequency would lower the upper bound on the stellar 
mass imposed by the radial instability limit, regard- 
less of the star's spin rate. 

Figure 3 shows clearly that observation of an SCO 
with a frequency > 1200 Hz usually reduces greatly 
the area of the radius-mass plane allowed for stars 
constructed using a given equation of state. For ex- 
ample, the area allowed for the UU equation of state 
if a 1400 Hz SCO is observed is only a small frac- 
tion of the area allowed by the requirement of stellar 
stability (see Fig. 3b). 

If the spin frequency of a star is known, the upper 
and lower bounds on its mass and radius imposed 
by observation of an SCO of a given frequency can 
be read off Figure 3 by looking for the intersections 
of the relevant bounding curves with the mass-radius 
curve for that spin frequency. For example, if the 
spin frequency is 600 Hz and the SCO frequency is 
1200 Hz, then a neutron star with the FPS equation of 
state must have a mass between 0.80 Mq and 1.82 Mq 
and a radius between 9.39 km and 12.26 km. 

Even if the spin frequency of a neutron star is un- 
known, one can still extract upper and lower bounds 
on the mass and radius of the star for a given equa- 
tion of state, using the highest observed frequency of 
an SCO from the source and the appropriate figure 
like Figure 3a or 3b. For example, if a 1200 Hz SCO 
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is observed, then a neutron star with the FPS equa- 
tion of state must have a mass between 0.64 M0 and 
2.12 and a radius between 9.28 km and 15.1 km. 

Figures 3a and 3b also show that for a given equa- 
tion of state, knowledge of any two of the mass, ra- 
dius, or spin frequency fixes the value of the third 
quantity, which may or may not be consistent with 
the kilohertz QPO frequency. Thus, if the spin fre- 
quency of a neutron star in an LMXB is known and 
the radius or mass can be determined by means other 
than observation of a kilohertz QPO (e.g., by measur- 
ing the emitting area during a thermonuclear X-ray 
burst or by measuring the mass dynamically), then 
observation of a kilohertz QPO will overdetermine 
the properties of the star, providing a check on the 
consistency of the mass and radius estimates. 

4. CONCLUSIONS 

Our results show that deviations from a first-order 
treatment of the effects of spin on the structure of neu- 
tron stars and on circular orbits around them are typ- 
ically significant for spin frequencies > 400 Hz. The 
Kerr spacetime is generally a poor approximation to 
the exterior spacetime of neutron stars spinning this 
fast or faster. 

Our results demonstrate that the upper bounds 
on the stiffness of neutron star matter implied by 
the high frequencies and coherences of the kilohertz 
QPOs are tightened significantly if the star is rotating 
rapidly. The constraints on the equation of state be- 
come much tighter if observations conclusively iden- 
tify a QPO frequency as the orbital frequency at the 
ISCO. For example, if the orbital frequency at the 
ISCO is 1100 Hz, the required mass is 2.0Mq for 
nonrotating stars and substantially more for rotat- 
ing stars. Such a high mass would rule out many of 
the softer equations of state and would imply that the 
three-nucleon interaction is strongly repulsive at high 
densities (see Pandharipande et al. 1998). 

As shown in § 3.3, even if the spin rate of the star 
is unknown, the area of the mass-radius plane allowed 
by observation of a high-frequency SCO can be quite 
small. If the spin rate is known, the range of radii 
allowed for a given equation of state is usually very 
small. 

Results of the kind presented in this work will be 
still more constraining if measurements of SCO fre- 
quencies can be combined with constraints on other 
quantities, such as the stellar compactness M/R (see 



Strohmayer 1997; Miller & Lamb 1998; Lamb et al. 
1998; Strohmayer et al. 1998) or the radius of the star 
(see Strohmayer et al. 1997b; Strohmayer et al. 1998). 

Additional observations of kilohertz QPO sources 
are extremely important, because these observations 
could provide very tight and robust constraints on 
the fundamental properties of neutron stars and on 
the equation of state of neutron star matter. 
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through the Compton Gamma-Ray Observatory Fel- 
lowship Program, by NASA grant NAS 5-2687. 

REFERENCES 



Akmal, A., Pandharipande, V. R., fc Ravenhall, D. G . 
1998, Phys. Rev. C. submitted ( |nucl-th/9804027|) 

Arnett, W.D., & Bowers, R. L. 1977, ApJ Suppl., 33, 
415 

Cook, G. B., Shapiro, S. L., & Teukolsky, S. A. 1992, 
ApJ, 398, 203 

. 1994a, ApJ, 422, 227 

. 1994b, ApJ, 424, 823 



Friedman, B., & Pandharipande, V. R. 1981, Nucl. 
Phys. A, 361, 501 

Glendcnning, N. K. 1992, Phys. Rev. D, 46, 1274 

Hartle, J.B., & Thornc, K. S. 1968, ApJ, 153, 807 

Heiselberg, H., Pethick, C.J., & Staubo, E. F. 1993, 
Phys. Rev. Lett., 70, 379 

Kluzniak, W., Michelson, P., & Wagoner, R. V. 1990, 
ApJ, 358, 538 

Komatsu, H., Eriguchi, Y., & Hachisu, I. 1989a, MN- 
RAS, 237, 355 

. 1989b, MNRAS, 239, 153 

Lagaris, I.E., & Pandharipande, V. R. 1981, Nucl. 
Phys., A359, 349 

Lamb, F.K., Miller, M. C, & Psaltis, D. 1998, in 
Nuclear Astrophysics, Proc. International Work- 
shop XXVI on Gross Properties of Nuclei and Nu- 
clear Excitations, ed. M. Buballa, N. Norenberg, 
J. Wambach, fc A. Wirzba (Darmstadt: GSI), 114 



( |astro-ph/980234q ) 



11 



Lightman, A. P., Press, W. H., Price, R. H., & Teukol- 
sky, S. A. 1973, Problem Book in Relativity and 
Gravitation (Princeton: Princeton University Press) 

Lorenz, C. P., Ravenhall, D. G., & Pethick, C. J. 1993, 
Phys. Rev. Lett. 70, 379 



Miler, M C , k ^,^imh, F K 1993, ApJ, 413, 143 



1996, ApJ, 470, 1033 

1998, ApJ Lett., in press (|astro-pli/9711325D 



Miller, M. C., Lamb, F. K., & Psaltis, D. 1998a, ApJ, 
in press 

. 1998b, in preparation 

Pandharipande, V. R., Akmal, A., & Ravenhall, D. G. 
1998, in Nuclear Astrophysics, Proc. Interna- 
tional Workshop XXVI on Gross Properties of Nu- 
clei and Nuclear Excitations, ed. M. Buballa, N. 
Norenberg, J. Wambach, & A. Wirzba (Darmstadt: 
GSI), 11 

Pandharipande, V. R., & Smith, R. A. 1975a, Nucl. 
Phys., A237, 507 

. 1975b, Phys. Letters, 59B, 15 

Pethick, C.J., & Ravenhall, D. G. 1995, Ann. Rev. 
Nucl. Par. Sci., 45, 429 

Schiavilla, R., Pandharipande, V. R., & Wiringa, 
R.B. 1986, NucL Phys., A449, 219 

Smith, D. A., Morgan, E.,H., & Bradt, H. 1997, ApJ, 
479, L137 

Strohmayer, T.E. 1997, talk presented at the 1997 
meeting of the High Energy Astrophysics Divi- 
sion of the American Astronomical Society in Estes 
Park, Colorado 

Strohmayer, T. E., Jahoda, K., Giles, A. B., & Lee, 
U. 1997a, ApJ, 486, 355 

Strohmayer, T.E., Zhang, W., & Swank, J.H. 1997b, 
ApJ, 487, L77 

Strohmayer, T., Zhang, W., Swank, J.H., Smale, A., 
Titarchuk, L., & Day, C. 1996, ApJ, 469, L9 

Strohmayer, T. E., Zhang, W., Swank, J.H., White, 
N.E., fc Lapidus, I. 1998, ApJ Lett., in press 
( |astro-ph/9803119| ) 

Swank, J. 1997, talk presented at the Symposium on 
The Active X-ray Sky: Results from BeppoSAX 
and Rossi-XTE, Accademia dei Lincei, Rome, Oc- 
tober 1997 



van der Klis, M. 1995, in X-Ray Binaries, ed. W. H. G. 
Lewin, J. van Paradijs, & E. P. J. van den Heuvel 
(Cambridge U.P.), 252 

. 1998, in The Many Faces of Neutron Stars, 



Proc. NATO ASI, Lipari, Italy, in press (astro- 
ph/97100iq ) 



Wijnands, R., Mendez, M., van der Klis, M., Psaltis, 
D., Kuulkers, E., & Lamb, F. K. 1998, ApJ, sub- 
mitted 

Wijnands, R.A.D., & van der Klis, M. 1997, ApJ, 
482, L65 

. 1998a, lAUC 6876 

. 1998b, submitted to Nature ( |astro-ph/9804216|) 



Wijnands, R. A. D., van der Klis, M., van Paradijs, 
J., Lewin, W. H. G., Lamb, F. K., Vaughan, B., & 
Kuulkers, E. 1997, ApJ, 479, L141 

Wiringa, R.B., Fiks, V., & Fabrocini, A. 1988, Phys. 
Rev. C, 38, 1010 

Zhang, W., Lapidus, I., Swank, J. H., White, N. E., 
& Titarchuk, L. 1997, lAU Circ. 6541 



This 2-column preprint was prepared with the AAS lAT^TjX 
macros v4.0. 



12 



Table 1: Mass-Radius Relations for the FPS Equation of State at Different Spin Rates 
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Table 2: Mass-Radius Relations for the UU Equation of State at Different Spin Rates 
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